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Abstract 

We investigate whether self-maintained vacuum traversible wormhole can exist 
described by stationary but nonstatic metric. We consider metric being the sum 
of static spherically symmetric one and a small nondiagonal component which 
describes rotation sufficiently slow to be taken into account in the linear approx- 
imation. We study semiclassical Einstein equations for this metric with vacuum 
expectation value of stress-energy of physical fields as the source. In suggestion 
that the static traversible wormhole solution exists we reveal possible azimuthal 
angle dependence of angular velocity of the rotation (angular velocity of the local 
inertial frame) that solves semiclassical Einstein equations. We find that in the 
macroscopic (in the Plank scale) wormhole case a rotational solution exists but 
only such that, first, angular velocity depends on radial coordinate only and, sec- 
ond, the wormhole connects the two asymptotically flat spacetimes rotating with 
angular velocities different in asymptotic regions. 
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1. Introduction. The possibility of existence of static spherically-symmetrical 
traversible wormhole as topology-nontrivial solution to the Einstein equations has been 
first studied by Morris and Thorne in 1988 Jl|. Since that time much activity has been 
developed in studying the wormhole subject (see, e.g., review by Visser [[|). Rather inter- 
esting is the possibility of existence of self-consistent wormhole solutions to semiclassical 
Einstein equations. Checking this possibility requires finding vacuum expectation value 
of the stress-energy tensor as functional of geometry and solving the Einstein equations 
with quantum backreaction, i.e. with such the induced stress-energy as a source. 

Recently some arguments in favour of this possibility has been given. In Refs. || f|, || 
the gravity induced vacuum stress-energy tensor in the wormhole background has been 
found to violate energy conditions just as it is required for this tensor itself be the 
source for such the wormhole metric |], ||. (We consider physical vacuum of spin 1 and 
1/2 massless fields in these papers). In Ref. 0] self-consistent spherically-symmetrical 
wormhole solution has been found numerically for the quantised scalar field vacuum 
playing the role of a source for gravitation. 

The problem of existence of self-consistent static wormhole is a particular case of the 
more fundamental problem of self-consistent solutions to semiclassical Einstein equations 
with vacuum expectations of stress-energy tensor of physical fields as a source. In Ref. 
II it has been found that such the problem linearised over metric perturbations off 
Minkowski spacetime gives solutions (apart from some unphysical ones) coinciding with 
those for classical gravity wave problem. If, however, topology is not Minkowski one as 
in the wormhole case at hand, new solutions can appear such as static wormhole itself. 
The natural next step may be the search for stationary but nonstatic topology nontrivial 
solution, namely rotating vacuum wormhole. In the case of slow rotation one simply adds 
small nondiagonal polar-angle-time component of metric to the self-maintained static 
wormhole metric: 

ds 2 = exp (2$)rft 2 - dp 2 - r 2 [d9 2 + sin 2 # (d(p 2 + 2hd<pdt)} (1) 

where h(p, 9) has the sense of angular velocity of the local inertial frame and will be 
called the angular velocity of rotation in what follows. The static metric $(p), r(p) is 
assumed to exist as a self-consistent solution of static semiclassical Einstein equations, 
and solution for h is to be found. One can take h arbitrarily small in order to limit 
oneself to the theory linearised in h. In practice, in order that the quadratic in h 
terms in the Riemann tensor could be disregarded, the space derivatives of h should 
be negligible as compared to the derivatives of the static part of metric in the scale of 
typical wormhole size. Note that in the linear approximation the only component of 
Einstein equations being new compared to the static case is the t<p one, as it follows 
from symmetry considerations. The role of the source is played there by the induced 
vacuum energy flow (angle-time stress-energy component). 

In the given note we show that the rotational solution exists at least in the case when 
the coefficient at the Weyl term in the induced stress-energy is large and if one can neglect 
other terms. As noted in Refs. [J, this corresponds to the wormhole of macroscopic 
size (that is, large in Plank units). The rotational solution which we argue should exist 
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is such that h depends on radial distance p only and has different asymptotic limits 
in the two asymptotically flat spacetimes connected by wormhole. This, in particular, 
means that these two spacetimes cannot be "glued" together in asymptotic region, so 
that this wormhole cannot be considered as that connecting the two regions of the same 
asymptotically flat spacetime. 

As for the general case when we do not assume the Weyl term to dominate, we 
find the two kinds of azimuthal angle 9 dependence of angular velocity h for which the 
radial and angular variables p, 9 are separated and the Einstein equation with quantum 
backreaction for h reduces to that for the function of purely p. One possibility is the 
above mentioned angle-independent angle velocity; another one is proportional to cos 9 
velocity. 

2. Classical rotation. By classical we mean rotation considered without taking 
into account corresponding quantum backreaction, i.e. induced vacuum energy flow. 
However, it is implied that the static wormhole problem is already solved (the static 
metric is found) with taking into account corresponding backreaction. We show (in the 
linearised in h theory used throughout the paper) that only solution for h not depending 
on 9 is physically acceptable such that it has different finite limits at p — > +oo and at 
p — > — oo. 

The t(f) Einstein equation of interest can be conveniently written using the tetrad 
components introduced like the following basic 1-forms u a = e^dx^ 0: 

oo = exp ($)dt, oo 1 = (d(j) + hdt)r sin 9, u 2 = dp, u 3 = rd9. (2) 

Taking the expressions for the Riemann tensor presented in Ref. || we find for the 
equation of interest: 

r . exp($) A Id , ^ A dh 




- $ )-^^ sin3 ^ = °- ( 3 ) 
sm z 9d9 89 

Remind that we neglect the induced vacuum energy flow T\ in this section. Separating 
the variables, h = f(p)Y(9), we find hypergeometric function for the angle dependence: 

y~F(a,3-a;2;^|^), (4) 

a = const. This function diverges at 9 = it unless a = —k, k = 0, 1,2, in the latter 
case it reduces to the Gegenbauer polynomial C^ 2 \cos9). Then the radial part satisfies 
the equation 

(r 4 exp (-$)/')' = k ( k + 3)r 2 exp (-$)/. (5) 

By asymptotical flatness $->0, ±oo and thus / ~ p~ k ~ 3 at large 

distances. Another solution, / ~ p k , should be discarded at k > 1 as unphysical one. 
Therefore if we choose some large L > 0, we shall have solutions at \p\ > L parametrised 
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by two constants C± (/ — > C±p~ k ~ 3 at p — > ±00), whereas in the intermediate region 
\p\ < L the Eq. (|5|) is regular in the wormhole geometry and has a solution parametrised 
by two constants C±, CV The four constants Ci, C2, C + , C_ are subject to four uniform 
equations which are matching conditions for / and it's derivative /' at p — +L and at 
p = —L. Requiring for this system to have nonzero solution, i.e. zero determinant, we 
get some constraint on the already known metric functions r(p), $(p). Therefore the 
set of possible solutions for f(p) at k > 1 has zero measure as compared to the set of 
possible static solutions r, $. At k = the Eq. (^) can be easily integrated and leads 
to physically admissible h not depending on 9 and being monotonic function of p which 
has finite but different limits at p — > ±00: h(—oo) 7^ h(+oo). Because of the latter 
circumstance the two asymptotically flat spacetimes connected by the wormhole channel 
cannot be "glued" together in asymptotic region (glueing at a shorter distance would 
spoil spherical symmetry of the background static solution described by two functions r, 
$) so we cannot derive from such the rotating wormhole the wormhole connecting the 
two distinct regions of the same asymptotically flat spacetime. 

3. Quantum backreaction and angle dependence of rotation. Now consider 
possible dependence of the angle velocity h on the azimuthal angle 9 that could solve 
the semiclassical Einstein equation (with backreaction). Here we show that the only two 
versions of the angle dependence for this equation to be solved by separation of variables 
are the following ones: h = f(p) or h = f(p) cos9. 

Evidently, the problem reduces to studying the angle dependence of T w for a given 
angle dependence of h. Given any physical field, we should solve equations of motion for 
this field in curved spacetime and sum vacuum contributions into stress-energy from all 
the eigenmodes. The resulting expression can be regularised by, e.g., covariant geodesic 
point separation and renormalised by subtracting the divergent parts known for physical 



fields PTJ|] . Choosing such separation in the radial direction we avoid discussing the 
renormalisation issue as far as the angle dependence is concerned. 

Let us consider general structure of the equations of motion and stress-energy for 
arbitrary field in the metric (HD and illustrate this by the case of massless fields of spin 
1 (electromagnetic) and 1/2 (neutrino). Most natural to display effect of rotation in 
the stationary axisymmetrical nonstatic metric is to use the complex Newman-Penrose 
formalism in analogy with that applied to Kerr metric || . In particular, up to the linear 
order in h, we can choose isotropic tetrad of real Z M , n M and complex m M , m*^ (asterics 
means complex conjugation) of which Z M , n M are tangential to some geodesies and m M 
(and thus m* M ) are orthogonal to / M , n M . Besides that, normalisation can be chosen such 
that / M n M = 1, m M m* = — 1. These are the distinctive properties of the Newman-Penrose 
tetrad which can be written with the help of the derivatives over directions as 



Pdp = exp (-2$)[d t + exp ($)d p - hd 4 
n^df, = -{d t -exp($)d p - , 
m"^ = (rV^y^de + i(sm9)- 1 i 



n^d^ = \[dt- exp ($)0 P - hd^ (6) 
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m*^ = (rV2)- l [d e - i(sm9)~%\. 

Calculation gives the following values for the standard 12 complex spin-connection coef- 
ficients k, u, a, A, e, g, /i, r, tt, a, /3, 7 ||: 

K = 1/ = 0, a = 2exp (-2$) A = -2e = - exp (— 2$)/i e sin 0, 

q = 2fi = -— exp (-$), - r = tt = ^=exp (-$)/i p rsin6», (7) 

1 cot % 1 z 

-a = /3= 7~y| ex P (-$) V" sin 9, 7 = exp ($) - -h e sin 

where subscript on h means corresponding derivative; prime on r, $ means derivative 
over p. The main feature of the covariant equations of motion for an arbitrary physical 
field is therefore occurence of h in the form h p sm8, hgsmO and hd^ there. In the 
diagrammatic language, /i-field-field vertex is combination of these expressions. Besides 
that, in the Newman- Penrose formalism operators acting on the angle variables appear 
in the form of spin raising and lowering operators 

C s = d e -^- + scot9, Ct = dg + ^- + scot9, (8) 
sin 6 sin a 

< s < so, So being spin of the field. These operators simplify in the basis of spin 
spherical harmonics proportional to the elements of rotation matrix D l sm ((j), 8, 0) C s 



and transform D l sm to D l s+lm and D l s _ lm , respectively. In this basis the expression 
for T10 turns out to be a value of spin weight s = ±1 (and m = 0), that is, combination 
of D\ Q (or D l _i 0) for different I. Quantum contribution to T w can be viewed as some 
loop diagram with external Tio- and h-legs. Since h can be expanded in the Legendre 
polynomials Pk ~ Dq , we can take _Dq as probe function for the angle dependence 
of h. Then the expressions h p sin 6 and hg sin 9 appearing in the vertex on /i-line are 
combinations of the values -D^ 1 and D^q 1 . This corresponds to the angular momenta 
k + 1 and fc — 1 flowing through the /i-line. It is less evident but shown at the end of 
this section that the vertex hd^ corresponds to the combination of angular momenta 
k + 1, k — 1, k — 3, ... . By conservation of angular momentum the T 10 also should 
be combination of Diq 1 ^ 271 , n = 0,1,2, ... . In particular, at k = 0, 1 the only term 
(D\ or D\ Q ) remains and Tio factorises into the functions of p and of 9. Moreover, since 
D\q 1 ~ Cl 12 just for k = 0, 1, the same ^-dependence also factors out in the LHS of 
Einstein equation. Therefore we conclude: h = f(p) or h = f(p)cos9 solves for the 
^-dependence of semiclassical Einstein equation. 

Finally, let us illustrate the above said by the examples of electromagnetic and neu- 
trino fields; for more detail on the Newman-Penrose description of these fields see Ref. 
0. Electromagnetic field is described by three complex functions / , fx, $2, for our choice 
of complex tetrad (Bp these are related to the electromagnetic field strength tensor F pu 
as follows: 



2/o = F w + hFe^ + exp ($)F pe + [F t4> + exp ($)F p4 



sin#' 
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2 f, = _ r 2 exp (-$)(F tp + hF p4> ) + ify^-, (9) 

sm u 

2/ 2 = -(Fte + hF 94> ) + exp + [F H - exp 

sm (7 

The eight real Maxwell equations can be recast into the following four complex ones: 
d-fi - Cifo = hd^fi, 

d + h + Ctf Q = -hd 4 ,f u (10) 
d+fo + CXP y^— Ctfi = i — 9 sin 6> - ih p exp ($) sin 6»/i - hd^fo, 

d-f 2 - 6Xp( 2 2$) £ /i = -i^^he sin + ifc, exp ($) sin 9f t + /ify/ 2 , 

where d± = exp ($)d p T- = w being energy. Each mode should be normalised so 

that it's full energy 

J T}r 2 exp ($)dp sin 0d0# = ^ sin 0<$d0 exp (-$)dp j/ 7o + /2V2 

, exp (2$) 



-/r/i-2/isin^m[(/ -/ 2 )7i] (H) 



be equal to the vacuum value uj/2 and then substituted into the expression for the 
stress-energy component studied, 

Tm L ^m T >\ = _ 2 5^Et£) Xm[(/o _ hYh] . (12) 

y r sm 9 v J r 6 

Analogously, massless fermion field is described by two complex values gi, g 2 obeying 
the field equations 

d -9i + 6XP ^^ C 1/2 g 2 = ^ sin9(g 1 h e - rh p g 2 ) + hd^, 
exp^$ £+ 29l - g + g 2 = % - sin 9{-rh p g x - g 2 h e ) + hd^g 2 . (13) 
Expression for the energy of each mode takes the form (on the field equations) 

J T/r 2 exp ($)dp sin 9d6d(f) = J sin 0d6>d</> exp {-®)dp ■ 4w(glg! + g 2 g 2 ) (14) 
while the stress-energy component of interest is 

Tl o = 6XP ^3 ^ {d*i^t/29i + {£>t/ 2 gi)*gi ~ 9*2^-1/292 ~ (Ci/ 2 g 2 )*g 2 - d e {g* 1 g 1 - g* 2 g 2 ) 

+ (r$' - r')(g* gi + g*g 2 ) + 2r exp (-<f>M(dt - hd^g, - g*(d t - hd^)g 2 ]} .(15) 



Solutions to the above equations of motion can be constructed iteratively. In zero 
order one takes 




/ d-RD 



— 1,771 



V 



1(1 + l)RD l 0m 
d+RDi 



(16) 



+ 1,771 



which are the well-known TE-modes for the electromagnetic field (TM-modes follow by 
multiplying this by i = ^/—l) and 



(0) 



Z 1 D[ 



1/2,771 



Z 2 D l _ 1/2tm 



(17) 



for the neutrino field. The R and Zi, Z 2 are some radial functions. Substituting Eqs. 
(16) and (|17D into the RHS of the equations of motion we find for the first 0(h) correction 



an expression of the type 

/M (1) 




/ . . . D{ 



1,771 



I k~ 

m 



J 

m 



1,771 



I k + 1 j ■ 
m — m 



+ 



n j 

■ -^ + 1,771 . 



m 



..Dl 

■ ■ ■ u §m 

■ ■ ly +l,m , 

I k J 

m — m 



and quite analogous one for the fermion field. The dots mean (real) factors which do 
not depend on 8, <fi and m. The linear order in h of interest comes from interplay in 
the bilinear T 10 between zero order ( |T6"D and first correction (0). By properties of 3j- 
symbols and rotation matrix elements D l sm summation over m just yields combination 
of the harmonics Diq P , p — 1, —1, —3, . . .. Important is that the last term in ( |T%| ) which 
stems from hd^ operator in the equations of motion is representable as combination of 
3j -symbols as 

1 k 3 \ = y A n ( l k + l ~ 2n 3 

m — m J n \m — m 



in 



(19) 



where A n does not depend on m. In Tio this and other terms enter multiplied by 
D l Qm (D 3 ^ l m )* and summed over m thus giving just combination of D^q 1 ^ 2 " 1 , n > 0. 

4. Macroscopic wormhole and radial dependence of rotation. Usually, if one 
does not assume existence of fundamental scales in the theory other than the Plank 
scale one expects the typical wormhole size be of the Plank scale too. However, a new 
scale can exist connected with coefficient of the Weyl term in the effective action. This 
coefficient is subject to renormalisation in both infrared (if massless fields are present 
in the theory) and ultraviolet regions. Possible large value of this coefficient can enable 
existence of the wormhole of macroscopic size. 

Here we argue that if the Weyl term coefficient is large and one can disregard other 
terms in the effective action then the conclusion concerning the existence of rotating 
wormholes resembles that for the case of classical rotation in Sect. 2. 
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The effective gravity Lagrangian density with taking into account the Weyl term can 
be written as proportional to R + (2p 2 )~ 1 C IJil/ \pC^ uXp Up to the full derivative, the 
Weyl term C ilv \ p C llvXp is equivalent to 2(R flu R pu — |-R 2 ). We calculate the latter up to 
the second order in h (required to get the first order in the equations of motion) using 
Riemann tensor given in Ref. in the tetrad components. Varying in h gives the desired 
^-component of the Einstein equations. Consider both versions of ^-dependence of h 
found in Sect. 3, h ~ 1 and h ~ cosd, and introduce new variable z via dz = exp (— $)dp 
and the function f = rexp (—<&). For h not depending on 9, h = f(p) = f(p(z)), the 
result reads 



f 4 /*exp (2$)), 



r 1 



+ 



10 



f 2 



1) + ~r zz r 



(subscript z means differentiation over z). For h = f(p) cos 9 we find 



f 4 /* exp (2$))* -4 exp (2$)f7 



+ 



10- 2 2 _ 
— r r H — r zz r 
3 z 3 



34 

y 



) f 2 f z 



r 4 



-r,,r 



r 2 f 2 



(20) 



(21) 



The infrared contribution to the coefficient /i~ 2 goes from the massless fields. In the 
considered case of fi~ 2 large in the Plank scale the typical wormhole size r 2 , is defined 
just by /i~ 2 as r 2 . = (3/x 2 ) -1 [§], ||. For example, in the vacuum of Ni spin 1 and N1/2 
spin 1/2 massless fields we have 



(3/i 2 



G 



1207T 



(4Ni + N 1/2 ) In 



/120tt A 2 



V g m 1 + N 1/2j 



(22) 



A being infrared cut off. 

Consider first Eq. ( |2"0"D which upon integrating both parts and denoting f z = g 
reduces to the second order one: 



~4 



/rexp (2$)r 4 s- + C = r 
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1) + -^f zz f 



f 2 g 



(23) 



where C = const. Asymptotical p — > ±00 form of this equation reads 



did 
dp p 4 dp 



P A 9 



c_ 



(24) 



The general solution is the sum of a particular one which behaves as g ~ p~ 4 + 0(p~ 6 ) 
at p — > ±00 and arbitrary combination of the two independent solutions to the uniform 
equation. Of the latter two one exponentially grows at p — >• +00 or at p — > —00 and 
should be omitted as unphysical solution while another one proportional to exp (—pp) 
(at p — > +00) or exp (pp) (at p — > —00) should be kept. Therefore, if we choose, as 
in Sect. 2, some large L > we shall have physically acceptable solutions at \p\ > L 
parametrised by three constants, one of which is C. Meanwhile, in the intermediate 
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region \p\ < L the equation ( PB| ) is regular in the wormhole geometry and has solution 
parametrised by maximal set of three constants, one of which is C . The overall set of 
five constants is subject to four uniform equations which are matching conditions for g 
and for it's derivative at p = +L and at p = —L. This defines all five constants up 
to an overall factor. Note that imposing additional condition gdz = (that is, 
h(—oo) = h(+oo)) is, generally speaking, contradictory since it would be condition not 
on a freely chosen constant, but on the already defined static metric $(p), r(p). 
Next consider Eq. fl2"Ip which has asymptotic form (at $ = 0, r = p) 

p 4±]_± p 4 f= ]_ p 4 (±]_± p *Y f (25) 

dp p 4 dp p? \dp p A dp ) 

Assuming this form at \p\ > L we find 



-f J J -fpV-/^ = C ± 1 P+# (26) 
dp p 4 dp p 4 



where Cj^, C+ 4 (CZ 1 , C+ 4 ) are some constants which parametrise the solution at p > L 
(at p < —L). To get physical solution we put C_ 1 = 0. Of the two solutions of uniform 
equation we discard exponentially growing one in each region p > L or p < —L and 
retain exponentially falling off another solution. Thus, in each region p>Lorp<—L 
the solution is specified by two constants. At the same time, the regular fourth order 
differential equation has solution parametrised by four constants at |p| < L. The overall 
number of constants is eight. These should ensure validity of eight matching conditions 
for /, /', /" and at p = ±L. The determinant of this uniform system should be zero. 
This imposes a constraint on the already known static metric $, r. Therefore the subset 
of rotating wormhole solutions with h = f(p) cos# should have zero measure w.r.t. the 
set of spherically symmetrical static wormhole solutions $(p), ?"(p). 

5. Conclusion. We have shown that if the coefficient at the Weyl term is large 
(infrared cut off is large) and one can discard other terms in the effective action then the 
rotation existing for any static wormhole background $(p), r(p) is that which proceeds 
with the angular velocity h not depending on the azimuthal angle 9 and having the 
different finite limits h(+oo) and h(— oo) in the asymptotic region p — > ±oo. We see that 
despite that the structure of the equations is drastically changed because of enhancing 
the maximal order of derivatives upon taking into account vacuum polarisation, the 
result practically does not differ from that for the classical case of Sect. 2. 

Note that the sign of infrared divergent coefficient p~ 2 at the Weyl tensor squared in 
the effective action is crucial for existence (or, rather, nonexistence) of more rotational 
solutions in our case. Were p~ 2 substituted by negative value, the equations above would 
have oscillating instead of monotonic exponential solutions, and we would not have to 
omit some of them as unphysical ones. Then the general solution of interest would be 
parametrised by more constants, and the set of such solutions would be larger. 

Also we can say that if we denote p? = x and extend the equations to arbitrary real 
x, the problem will be singular at x = 0. 
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In the case if the infrared logarithm is not large we do not have macroscopic vacuum 
wormhole, and the following interesting question arises: whether microscopic wormhole 
can rotate so that it would have the macroscopic "tail" of rotation (when h falls off in 
power law in asymptotic region). Answering this question implies rather complicated 
problem of calculation and analysis of the terms in stress-energy other than the Weyl 
term. 

This work was supported in part by the President Council for Grants through grant 
No. 96-15-96317. 
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